NASA Technical Memorandum 108492 





ANSYS Duplicate Finite-Element Checker Routine 


R. Ortega 


(NASA-TM-108492) ANSYS DUPLICATE 
FINITE-ELEMENT CHECKER ROUTINE 
(NASA. Marshall Space Flight 
Center) 20 p 


G3/6^ 0055A63 


N95-2982A 
Unci as 


May 1995 




NASA Technical Memorandum 108492 


ANSYS Duplicate Finite-Elem ent Checker Routine 

R. Ortega 

Marshall Space Flight Center • MSFC, Alabama 


National Aeronautics and Space Administration 
Marshall Space Flight Center* MSFC, Alabama 35812 


May 1995 




TABLE OF CONTENTS 


Page 


I. 

INTRODUCTION 

1 

n. 

THEORY 

1 

m. 

ROUTINE DESCRIPTION 

0 

IV. 

ATD HPOTP RESULTS ... 

Q 

V. 

CONCLUSIONS 

0 


REFERENCES 

APPENDIX - Duplicate Elements Checker ANSYSSOa Routine 


t;C£D|flQvFAQC BLAfJK NOT ris.MED 
page 4I mrnimmj 


iii 


LIST OF ILLUSTRATIONS 

Figure 

1. Common face of adjoining elements 

2. Planar view of points A, B, C, and D 

3. Examples of floating elements within a mesh 

4. Examples of fully duplicated elements 

5. Examples of intersecting elements with common faces .... 


Page 


5 

5 

6 

7 

8 


IV 



TECHNICAL MEMORANDUM 


ANSYS DUPLICATE FINITE-ELEMENT CHECKER ROUTINE 

I. INTRODUCTION 


This report documents an ANSYS routine to check for duplicated elements within the volume 
ot a three-dimensional (3-D) finite-element mesh.i Recent finite-element submodeling efforts of the 
space shuttle mam engine alternate turbopump development (ATD) high pressure oxidizer turbo- 
p^ump (HPOTP) uncovered duplicate elements in a certain location of the global finite-element model 

^ ATD HPOTP which consists of a total of 

54,296 elements distributed over 14 superelements. An automated method was needed to check the 
large number of elements involved. 2 A routine was then developed to check 3-D finite-element 

meshes for the existence of unconnected elements, elements within elements, and elements across 
other elements. 


II. THEORY 


Figure 1 shows the common face of adjoining elements (face BEGF) and the tetrahedral 
volumes formed by Joining the corners of this common face to the respective element centroids 

^ intersection of the line passing through A and D and face 

provide the coordinates of point C 

(AC, YC, ZC). Tnis IS done mathematically by defining the following cross products: 


RA = ABxAD , 


RB = AExAD , 
^ = bBxbP . 


Vectors RA, RB, and RC represent the normal vectors of the three planes defined above. 
Using these normal vectors, the equations of the three planes can be written as a system of three 
linear equations; ^ 


aaXC+abYC+ac7£ = PA , 
baXC+bbYC+bcZC = PB . 
caXC + cbYC + ccZC = PC . 


Cram^^ rule, found in any algebra book, can be used to solve the above system providing coordi- 
nates AC, YC, and ZC of point C.^ 



Figure 2 shows points A, B, C, and D. The line joining points B and C falls within face BEGF 
(shown in fig 1). For adjoining elements not to have any part of their volume crossing each other, it 
Semonslraced thil points A and D (element centroids) fall on oppostte stdes of the hne 

joining B and C. 

This can be done by comparing the normal vectors from the planes obtained from ^eas 
and DBG. The dot product of two unit vectors provides the cosine of the be^ween^^^^^ 

If the value is negative, then the angle must be greater than 90 

dot product of the normal vectors from the planes obtained from areas ABC and DBC m^t be less 
than zero for the centroids A and D to fall on opposite sides of the line joining B and C. Thus, the 
following inequality must be satisfied: 

(ABxAC) (DBxDe)<0 . 

Otherwise, the elements are deemed to have a volume intersection with each other. There 
are two key assumptions that must be noted. First, for elements to be adjoining each other, they 
must have a face in common, which is accomplished by sharing at least three of the four comer 
nodes. Second, the face of adjoining elements is considered to be flat or with a very minor degree 

warping. 


III. ROUTINE DESCRIPTION 


To implement the check of the above develop inequality an ^S^^Oyouto ^ 
ten and is enclosed in the appendix. In order to mn this routine, a data base file (.db file) of the model 
in ANSYSSOa format must be available. The routine is called checker3d3^at and runs in a c mo . 
The data base file of the job that is to be checked should be copied into a file named duplicate.db and 
placed in a folder together with checker3d3.dat. 

The ANSYSSOa routine first resumes the model from the data base file and counts the nuin- 
ber of elements (numb). An array is then set with a size of numb by six. The number six represents 
the maximum number of sides that a brick element is allowed to have. A do loop is st^ed to check 
evereleTem on every side. Because of this, the routine takes a very large amount of central 
proSssing unit (CPU) time to complete its task. The routine was onginally developed as an 
SsYSSOa macro with the intent of using ANSYS select logic. The assumption was that using ^ 
already developed select logic would reduce the time necessary for developing a routine to conduct 
L element chLks. The amount of CPU runtime savings to be obtained from developing a checker 
FORTRAN routine from scratch (including select logic) is not known, however, it is suspected that 

the savings could be quite large. 

Three element irregularities are checked by the ANSYSSOa routing The first is when 
an element without any other elements attached to it at a common face. This means that the element 
could be floating within the mesh attached to other element nodes. Figure 3 shows examples of this 
condition Duplicate elements are shown in gray in the figure. Unattached floating elements are 
reported in a file called single.elem. The second element irregularity is when two elements are 
duplicated having exactly the same nodes. Fully duplicated elements numbers, shown m figure 4, 
are^given as pairs^in nontriv.elem. Fully duplicated elements are those that contain e^ctly the same 
nodes. Finally, the third condition checked is when partial or full volumes of element w^^ 

cross other elements and have a common face. The third condition is shown in figure 5. The eleme ts 
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numbers are reported in double.elem in pairs. Elements shown in double.elem should be visually 
checked by pairs, because element warping could result in the reporting of legitimate elements as 
intersecting each other. This is due to the fact that a warped face could wrap around a centroid 
location, resulting in a positive value for the above inequality. 


IV. ATD HPOTP RESULTS 


The results of the check for duplicate elements on the ATD HPOTP global model are shown 
in table 1. The table shows the superelement id, number of elements, the number of duplicate ele- 
ments reported by routine checker3d3.dat, and the final number of duplicate elements after the final 
visual check of the reported elements. The difference between the reported and final number of dupli- 
cate elements represents those elements whose degree of warping was such that the check inequal- 
ity was not satisfied, but are acceptable to the ANSYS code. The last column of table 1 reports the 
CPU time required to run the check on a Silicon Graphics ONYX workstation with two 100-mHz 
processors. 


V. CONCLUSIONS 


An ANSYS5.0 routine to check for duplicated elements within the volume of a 3-D finite- 
element mesh was developed. This routine was used to identify ATD HPOTP finite-element model 
duplicate elements. It is of importance to note that the number of duplicate elements is small com- 
pared to the overall model and ended up having little impact on the final results of the ATD HPOTP 
model. However, these duplicate elements could have occurred in a stress concentration region pro- 
viding erroneous stress results. Also, the mere existence of such elements can complicate 
submodeling efforts in regions where they are contained, as it did occur on a submodel effort for the 
ATD HPOTP weld No. 2. 

The ANSYS finite-element code assumes that duplicate elements are modeled on purpose. 
There are instances where duplicate elements are desirable in a mesh. However, unwanted duplicate 
elements might appear when a complicated 3-D geometry mesh is modified by hand. Therefore, it is 
recommended that models be checked when there is knowledge or suspicion of extensive use of 
manual node and element definitions within a model. Finally, considering the large amount of CPU 
time needed to run this check, it is suggested that the duplicate element check run be done on a copy 
of a data base file on a separate folder in the background while work on the model continues as if no 
duplicate elements exist. 
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Table 1. Duplicate elements check of ATD HPOTP. 


Superelement ID 
Number 

Number of 
Elements 

Number of 
Duplicate 
Elements 
Reported 

Final Number 
of Duplicate 
Elements 

CPU Time 
(h) 

SEOl 

3,456 

0 

0 

31.05 

SE 02 

7,132 

32 

28 

79.29 

SE 03 

4,007 

5 

0 

44.26 

SE 04 

3,168 

0 

0 

30.86 

SE 05 

3,484 

0 

0 

33.16 

SE 06 

2,924 

0 

0 

19.47 

SE 07 

3,254 

0 

0 

20.28 

SE 08 

7,227 

24 

18 

75.39 

SE 09 

5,972 

0 

0 

54.08 

SE 10 

2,267 

2 

0 

24.21 

SE 11 

4,420 

17 

15 

26.27 

SE 12 

5,414 

0 

0 

84.04 

SE 13 

687 

0 

0 

7.57 

SE 14 

884 

0 

0 

3.25 






















































Figure 3. Examples of floating elements within a mesh. 



Figure 4. Examples of fully duplicated elements 
(Two elements with identical nodes). 
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Figure 5. Examples of intersecting elements with common faces. 
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/filename.duplicate 

resume 

/prep7 

arr(l,l)= 

yes 

esel.none 

cra,crset,elem 

esel.all 

nsel,all 

cm, dummy ,elem 

*get,numb,elem,0,count 

*get,numbb,elem,0,num,max 

*dim,arr„numbb,6 

*do,in,l,numb,l 

*get,en,elem,0,num,min 

esel,s,elem„en 

nsle 

♦get,nn,node,0,count 

*if,nn,le,4,then 

nef=nelem(en,5) 

*if,nef,le,0,then 

cmsel.s.dummy 

esel,u,elem„en 

cm,dummy,elem 

♦cycle 

♦endif 

♦endif 

*get,xa,elem,en,cent,x 

*get,ya,elem,en,cent,y 

*get,za,elem,en,cent,z 

imark=0 

*do,jn, 1,6,1 

*if,jn,eq,l,then 

b=nelem(en,l) 

e=nelem(en,4) 

f=nelem(en,2) 

♦endif 

*if,jn,eq,2,then 

b=nelem(en,2) 

e=nelem(en,6) 

f=nelem(en,l) 

♦endif 

*if,jn,eq,3,then 

b=nelem(en,3) 

e=nelem(en,7) 

f=nelem(en,2) 

♦endif 

*if,jn,eq,4,then 

b=nelem(en,4) 

e=nelem(en,8) 


f=nelem(en,3) 

♦endif 

*if,jn,eq,5,then 

b=nelem(en,l) 

e=nelem(en,5) 

f=nelem(en,4) 

♦endif 

*if,jn,eq,6,then 

b=nelem(en,6) 

e=nelem(en,7) 

f=nelem(en,5) 

♦endif 

nsel,s,node„b 

nsel,a,node„e 

nsel,a,node„f 

cm,stud,node 

*get,nnb,node,0,count 

*if,nnb,lt,3,then 

imark=imark+l 

*if,imark,eq,6,then 

esel,s,elem„en 

ewrite,single,elem„ 1 

cmsel,a,crset 

cm,crset,elem 

♦endif 

♦cycle 

♦endif 

esel,none 

cm,elgr,elem 

nsel,s,node„b 

esln 

♦get,elb,elem,0,count 

cm,elst,elem 

♦do,jk,l,elb,l 

♦get,est,elem,0,num,min 

esel,s,elem„est 

nsle 

cmsel,r,stud 

*get,chk,node,0,count 

♦if,chk,eq,nnb,then 

cmsel,a,elgr 

esel,u,elem„en 

cm,elgr,elem 

♦endif 

cmsel,s,elst 

esel,u,elem„est 

cm,elst,elem 

♦enddo 

cmsel,s,elgr 

*get,elc,elem,0,count 
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*if,elc,le,0,ihen 

imark=imark+l 

*if,imark,eq,6,then 

esel,s,elem„en 

ewrite,single,elem„l 

cmsel.a.crsel 

cm.crset.elem 

*endif 

♦cycle 

*endif 

*if,elc,eq,l,then 
*get,ena,elem,0,num,niin 
ittt=0 
*do,jp, 1,6,1 

*if,arr{ena,jp),eq,en,then 

ittt=l 

*endif 

*enddo 

*if,ittt,eq,l,then 

*cycle 

*endif 

*endif 

*do,jl,l,elc,l 

*get,ena,elem ,0,num,min 

esel,s,elem„ena 

esel,a,elem„en 

*get,xd,elem,ena,cent,x 

*get,yd,elem,ena,cent,y 

*get,zd,elem,ena,cent,z 

xb=nx(b) 

yb=ny(b) 

zb=nz(b) 

xe=nx(e) 

ye=ny(e) 

ze=nz(e) 

xf=nx(f) 

yf=ny(0 

zf=nz(0 

al=xb-xa 

a2=yb-ya 

a3=zb-za 

ael=xe-xa 

ae2=ye-ya 

ae3=ze-za 

adl=xd-xa 

ad2=yd-ya 

ad3=zd-za 

bel=xe-xb 

be2=ye-yb 

be3=ze-zb 


bfl=xf-xb 

bf2=yf-yb 

bf3=zf-zb 

aa=a2*ad3-a3*ad2 

ab=a3*adl-al*ad3 

ac=al*ad2-a2*adl 

ba=ae2*ad3*ae3*ad2 

bb=ae3 *ad 1 -ae 1 *ad3 

bc=ael *ad2-ae2*adl 

ca=be2*bf3-be3*bf2 

cb=be3*bfl-bel*bf3 

cc=bel*bf2-be2*bfl 

pa=aa*xa+ab*ya+ac*za 

pb=ba*xa+bb*ya+bc*za 

pc=ca*xb+cb*yb+cc*zb 

mdl=aa*bb*cc-aa*cb*bc-ba*ab*cc 

md2=ba*cb*ac+ca*ab*bc-ca*bb*ac 

md=mdl+md2 

mxl=:pa*bb*cc-pa*cb*bc-pb*ab*cc 

mx2=pb*cb*ac+pc*ab*bc-pc*bb*ac 

mx=mxl+mx2 

myl=aa*pb*cc-aa*pc*bc-ba*pa*cc 

my2=ba*pc*ac+ca*pa*bc-ca*pb*ac 

my=myl+my2 

mzl=aa*bb*pc-aa*cb*pb-ba*ab*pc 

mz2=ba*cb*pa+ca*ab*pb-ca*bb*pa 

mz=mzl+mz2 

*if,md,eq,0,then 

cmsel,u,crset 

*get,nca,elem,0,count 

*if,nca,gt,0,then 

esel,s,elem„en 

esel,a,elem„ena 

ewrite,nontriv,elem„ 1 

cmsel,a,crset 

cm,crset,elem 

*endif 

cmsel,s,elgr 

esel,u,elem„ena 

cm,elgr,elem 

♦cycle 

♦endif 

xc=mx/md 

yc=my/md 

zc=mz/md 

bl=xc-xa 

b2=yc-ya 

b3=zc-za 

cl=xb-xd 

c2=yb-yd 
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c3=zb-zd 

dl=xc-xd 

d2=yc-yd 

d3=zc-zd 

0ne=(a2*b3-a3*b2)*(c2*d3-c3*d2) 

two=(a3*bl-al*b3)*(c3*dl-cl*d3) 

three=(al *b2-a2*b 1 )*(c 1 *d2-c2*d 1 ) 

fin=one+two+three 

*if,fin,ge,0,then 

cmsel.u.crset 

*get,nca,elem,0,count 

*if,nca,gt,0,then 

esel,s,elem„en 

esel,a,elem„ena 

ewrite,double,elem „ 1 

cmsel,a,crset 

cm.crset.elem 

*endif 

♦else 

arr(en,jn)=ena 

♦endif 

cmsel,s,elgr 

esel,u,elem„ena 

cm.elgr.elem 

♦enddo 

♦enddo 

cmsel,s,dummy 

esel,u,elem„en 

cm, dummy, elem 

♦enddo 

flni 
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